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ABSTRACT 

The Fe Ka emission line is the most ubiquitous feature in the X-ray spectra of active galactic nuclei (AON), 
but the origin of its narrow core remains uncertain. Here, we investigate the connection between the sizes of 
the Fe Ka core emission regions and the measured sizes of the dusty tori in 13 local Type 1 AGN. The observed 
FeKa emission radii (i?Fe) are determined from spectrally resolved line widths in X-ray grating spectra, and 
the dust sublimation radii (i?dust) are measured either from opticaPnear-infrared reverberation time lags or 
from resolved near-infrared interferometric data. This direct comparison shows, on an object-by-object basis, 
that the dust sublimation radius forms an outer envelope to the bulk of the Fe Ka emission. i?Fe matches i?dust 
well in the AGN with the best constrained line widths currently. In a significant fraction of objects without 
a clear narrow line core, i?Fe is similar to, or smaller than the radius of the optical broad line region. These 
facts place important constraints on the torus geometries for our sample. Extended tori in which the solid 
angle of fluorescing gas peaks at well beyond the dust sublimation radius can be ruled out. We also test for 
luminosity scalings of i?Fe, finding that Eddington ratio is not a prime driver in determining the line location in 
our sample. We discuss in detail potential caveats due to data analysis and instrumental limitations, simplistic 
line modeling, uncertain black hole masses, as well as sample selection, showing that none of these is likely 
to bias our core result. The calorimeter on board Astro-H will soon vastly increase the parameter space over 
which line measurements can be made, overcoming many of these limitations. 

Subject headings: galaxies: active — X-rays: galaxies — Infrared: galaxies 


1. INTRODUCTION 


The Fe Ka line is the most prominent emission feature ob¬ 
served in the spectra of accreting systems. It arises as a re¬ 
sult of fluorescence in predominantly cold gas, and has been 
observed in X-ray binaries as well as active galactic nuclei 
(AGN). The broad component has been modeled as origi¬ 
nating in th e inner accretion disk as a result of relativistic 
broadening (iFabian et al.lll989]i . Despite many years of study, 
though, the origin of the narrow core of the line in AGN is 
still unclear. Previous works have found no strong correlation 
of the Fe line width with those seen in the broad line region 
(BLR), and suggest that there are likely multiple sites of ori¬ 
gin, including the dusty torus, the BLR, as well as the outer 
accretion disk j Yaqoob & Padmanabhanll2004t iNandral120061 : 

iShii etal.ll2O10i). 


All modern AGN X-ray radiative transfer models include 
fluorescence emission computed self-consistently for a 
variety of geometries and obscuring colu mn densities in the 
accre tio n disk as wel l as the to rus (e.g. iGeorge & Fabi^ 
19911 iNandraetaT 


g007l: iMurnhv & YaaoobI 12009^ 
1201 Ih . and recent high quality X- 
ray spectra for many AGN - obscured as well as unobscured 
bear out the correspondence between transmission . 


Brightman & Nandra 


Compton scattering and fluorescence 
Brenneman et al.l 20141 [Arevalo^e^al 


2014t iBaueret al. 12014 Gandhi et al 


( e.g. r^aoobllini 

jiour iGandhi et alJ 

l2015h . Luminosity 


and covering factor scalings between the infrared and X-ray 
regimes may also support a close connection between the 
line-emitting gas and dis t ribution of warm dust jGandhi et alj 
[200l iRicci et al.l IMl iToba et alJ |20ll, but any spatial 
correspondence between the FeKa emission zone and the 
torus still remains to be tested in detail. 


The main issue hindering progress in pinning down the ori¬ 
gin of the line core is lack of high quality data. In the X-rays, 
the broad and narrow components of the line must be disen¬ 
tangled, and the narrow core is only resolvable using grat¬ 
ing spectra in the brightest AGN at present. In the infrared, 
the dust emission zones span sub-pc to pc size scales (e.g. 
Kishimoto et al.ll201 lat iTristram et al.ll20^ iBurtscher et alJ 
2013h . which require interferometric techniques in order to 
resolve directly. Infrared reverberation mapping is a grow¬ 
ing held a nd provides an alternate means to infer dust emitter 
sizes (e.g. ISuganuma et alJl2()()^lKoshida et al.ll2014l) . 

There have been recent advances on all these fronts, with 
a growing number of sources with commensurate data now 
available. Here we present a comparison of the directly re¬ 
solved inner edge of the dusty tori (using either interferometry 
or reverberation mapping) with the sizes of the narrow Fe Ka 
line core regions for a sample of 13 local AGN. We inves¬ 
tigate the detailed connection between the two in individual 
objects, discuss the limitations in the current data, and make 
comparisons to other recent works. We also test whether the 
AGN Eddington accretion rates play a role in determining the 
Fe Ka location. Our study, albeit on a small sample, sets the 
stage for much larger studies w hich will be possible o nce the 
calorimeter on board Astro-H (iTakahashi et al.ll2014ft begins 
operation in 2016. 

2. SAMPLE 

2.1. X-ray grating spectra 

Our star ting samp l e is th e ^High Energy Grating (HEG) 
sample of IShu et al] (1201 Ol) . the largest sample to date of 
high spectral resolution data covering the ^6.4keV energy 
regime. This includes 36 unique galaxies below z = 0.3. In 
27 of these, a measurement of the velocity full widths at half 
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maximum (ufwhm) of the narrow line core was reported by 
the authors. This includes positive FWHM measurements, as 
well as unresolved line upper limits. The mean ufwhm for 
this sample was found to be 2060± 230kms“^. The HEG 
spectral resolution is 0.012 A, corresponding to a velocity res¬ 
olution of « 1860km s“^ at 6.4keV, and is accounted for in 
the line modeling. 

We have chosen measurements that are likely to best rep¬ 
resent the narrow line core. In this respect, note that several 
new observa tions for NGC 40 51 have become available since 
the work of iShu et al.l (1201^. We use the res ults discussed 
by iLobban et al.l (1201 11) and IShu et al.l (120121) . The former 
work found that the presence of unresolved line component at 
6.4keV significantly improved the ht of the line complex in 
the time-averaged data from 12 observations totaling ^ 300 ks 
of exposure. Since this component is the one most likely as¬ 
sociated with distant material, we treat the core as being un¬ 
resolved. This leads to an upper limit on ufwhm and a corre¬ 
sponding lower limit on i?Fe- 

IC4329A is also treated as upper limit on wfwhm- A1- 
though a broad compon ent to the line is detecte d by 
I Shu et al] (1201 Oi see also iMcKernan & YaaoobI 120041) with 
wfwhm = 1 8830 + 9020 ° km s“^, the line peak energy in this 
case is found to be E'peait = 6.305+Q|!)ggkeV. An alternative 
fit with an unresolved line by the authors yielded a more plau¬ 
sible value of TJpeak = 6.399t°;°05 keV closer to the expected 
neutral line energy, and the authors stress that these are more 
reliable measurements of the true narrow core. 

For measurement of the Fe Ka emission radii (i?Fe), we as¬ 
sume virial motion for the emitting clouds. Then, 


Rfb 


GMbh 


( 1 ) 


where Mbh is the black hole mass, and u = •\/3/2vfwhm- 
The correction factor of '/3/2 arises under the assumption of 
an isotropic velocity distribution and accounting for the fact 
that the line-of-sight velocity dispersio n is half of the FWHM 
(cf. lNetzeil[T990t Peterson et al.ll2004ft . Uncertainties related 
to unknown geometric projection factors are discussed below 
and in the Appendix. Where possible, Mbh values based 
upon reverberation mapping measurements are used, mostly 
based upon the H/3 time lag. This was the case for all sources 
except IRAS 13349-)-2438 which is derived from constraints 
on modeling its spectral energy distribution. In this case, we 
assumed a large uncertainty equal to the mass measurement it¬ 
self in order to account for the lack of reverberation mapping. 
In addition, for NGC4151, we use a recent measurement ob¬ 
tained from dynamical modeling. References are given in Ta¬ 
ble 1. Uncertainties on i?Fe were determined by joint Monte 
Carlo sampling of Mbh and v in log space. 

One systematic uncertainty in converting velocities to sizes 
comes from the unknown geometric projection, which may 
result in underestimation of the true space velocities, leading 
to overestimated sizes. This is a well-known effect in BLR re¬ 
verberation mapping and single- epoch estimates of M rh, and 
is captured in the /-factor (e.g. iPeterson et al.ll20d^ . If we 
posit that conservation of angular momentum of the accret¬ 
ing matter leads to a flattened geometry, then we can expect 
that the projection effects are very similar for the X-ray emit¬ 
ting region and the BLR. In this case, size measurements of 
these regions will be affected in a similar way for any individ¬ 
ual object, so comparisons between these quantities will not 


be affected. On the sample level, the distribution of geomet¬ 
ric inclinations will lead to a widening of any size-luminosity 
relation above the intrinsic dispersion by a factor of the or¬ 
der unity. Given that our sample covers four decades in lu¬ 
minosity and two decades in sizes, we can expect that the 
essence of the relation will be preserved even for unknown 
projection effects in individual sources. In the Appendix, we 
will discuss alternate size comparisons using reverberation- 
independent Mbh values - a test which avoids the above as¬ 
sumption. 

2.2. Infrared and optical data 

Near-IR (NIR) observations probe the emission region of 
the sublimation zone in the torus at temperatures of ^ 1500 K. 
The size of this region can be inferred either from reverbera¬ 
tion mapping or directly from interferometry. Although both 
types of measurements are qualitatively similar, the small 
quantitative differences between bot h are probably related 
to the detailed dust distribution (e.g. iKishimoto et al.l 120091 
Honig & Kishimot^ 1201 It iKoshida et al.l l2014t iHonig et alJ 
2014ft . Here we collect and use both types of data from the lit¬ 
erature (collectively referred to as i?dust)- We note that i?dust 
measurements are available mostly for Type 1 AGN, because 
in Type 2 AGN the innermost hot dust is not easily visible. 

The dust reverberation mapping radii (i?dust,rev) are the 
result of a long-ter m monitoring c ampa ign with the MAG¬ 
NUM telescope by IKoshida et al.l (12014ft and have been in¬ 
ferred from V- and AT-band light curve time lags (rdust) as 
f?dust,rev = crdust/(l + z) where c denotes the speed of light 
and the (1 -f z) factor corrects for cosmological time dila¬ 
tion. The interferometric sizes (i?dust,intf) are based on Keck 
Interferometer data and represent ra dii of a thin-ring model 
(IKishimoto et al.ll2009Ll20lIil2Ml . 

Of the 27 so urces with vfwhm from the sample of 
IShu et all (120101) . 7 have i?dust.intf measurements, and 10 
have published Tdust values. Our final sample of objects with 
measured values of i?Fe as well as i?dust comprises 13 unique 
sources. These are listed in Table 1. 

For these objects, we also compute the radii of their opti¬ 
cal BLRs (i?BLR) for comparison to i?Fe- Measurements of 
the reverberation time lag of the H/3 emission line (th/s) for 
12 AGN of our sampl e - i.e. all except IRAS 13349-)-2438 - 
have been tabulated in lBentz et al.l (12009ft . which we use here. 
Then, Rblr = cth/3/(1 + z). 

We follow this by investigating any possible relation be¬ 
tween the locations of emitting regions and the contin¬ 
uum emission. As a proxy of the latter, we consider the 
monochromatic continuum lumin osities at 5500 A 
whose values a r e gath ered from IKis himoto et al.l (1201 lal) . 
ISuganuma et al.l (12006ft and iBentz et al.l (120091) ' TTiese are 
based upon fitting of their spectral energy distributions and 
corrected for starlight contamination. A mean uncertainty of 
0.1 dex is assumed for ^ 5500 - These are also listed in Table 1. 

3. RESULTS 

Fig. [U shows the comparison between i?Fe and i?dust- For 
all objects, the measurements or limits on i?dust,intf lie a lit¬ 
tle above i?dust,rev For sources with measurements of both 
radii, i?dust,intt is larger than i?dust,rev by an average factor 
of 2.3 (± 0.3 mean standard deviation). Reverberation is sen¬ 
sitive to the fastest dust response with changing incident radi¬ 
ation, whereas the emission probed in interferometry is more 
sensitive to the average emitting surfa ce area which is likely 
to peak at somewhat larger radii (e.g. IKishimoto et'ani2009l : 
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jHonig & Kishimot^ 1201 ll : iKoshida et alJ l2014t iHonig et aH 
[2014lh In any case, the trends discussed below are similar 
for both measures of i?dust- 

For 9 of the 13 AGN, the measurements or lim¬ 
its on i?Fe are fully consistent with i?dust- These are 
NGC3783, NGC4151, NGC4593, NGC5548, Mrk509, 
3C273, IRAS 13349H-2438, NGC4051 and IC4329A. The 
remaining 4 sources (NGC3516, NGC7469, Mrk590 and 
Fairall9) have i?Fe values significantly smaller than i?dust- 
The ratio of i?dust / Rpe in these 4 ranges over 2.5 
(NGC 3516) to 75 (Fairall9). 

The figure also plots i?BLR on the right-hand axis and re¬ 
veals some interesting results for individual sources. Whereas 
the median value of i?BLR is 8 times smaller than i?Fe for 
the full sample, there are 6 objects (NGC 4051, NGC 4593, 
NGC 7469, Mrk509, Fairall 9 and 3C 273) for which the i?Fe 
measurement or limit is entirely consistent with i?BLR- For 
NGC 3516, i?Fe lies below i?dust and the uncertainty esti¬ 
mates on i?BLR and i?Fe do not overlap, implying an Fe Ka 
origin in an intermediate zone between the BLR and the torus. 
On the other hand, for Mrk 590, i?Fe is significantly smaller 
than even i?BLR- Finally, for Mrk509 and Mrk590, i?BLR 
lies within a factor of 1.5 of i?dust, implying a close proxim¬ 
ity of the BLR clouds with the inner extent of the torus. 

Examining the overall distribution of sources, the most im¬ 
portant feature is the absence of sources significantly above 
the line of equality. Whereas there is significant scatter of 
sources to small values of i?Fe, for no source is i?Fe much 
greater than i?dust- The two lower limits are also consistent 
with this line. We discuss these results at length in the next 
section. 

We next test whether the location of the Fe Ka line is driven 
by fundamental luminosity scalings. We first plot i?Fe as a 
function of L 5500 in Fig.|2] This shows a close correspon¬ 
dence with Fig. [U in all aspects, with an absence of sources 
on size scales asso ciated with tho s e above the torus size- 
luminosity relation (lBarvainislll987l : [Kishimoto et al.ll201 lal) . 
and a strong scatter of sources below. We tested for any rela¬ 
tion between the Rpe and L 5500 with Monte Carlo resampling 
and computing the distribution of Spearman rank coefficients 
(p) for the randomized ensembles. We find p = 0.45lg'25 
(90%) and insignificant p-values over the entire ensemble, 
suggestive of only a weak positive correlation. Any stronger 
trend is hidden by the large scatter in the current (small) sam¬ 
ple. 

In Fig.|3l we remove black hole mass scaling from both 
axes by plotting i?Fe in units of the Gravitational ra¬ 
dius (Rg =GMbh/c^) as a function of the Eddington ratio 
(TBoi/TEdd)- The latter quantity is taken as a proxy of 
the (specific) accretion rate. Bolometric luminosities (Lboi) 
are approxima te estimates based upon a correction factor 
of 6 to L 5500 (IScott et al.ll2004t [Richards et al.ll2006h . We 
cross-checked these LboI values with those reported from 
the bro adband modeling carried out by IVasudevan & Eabi^ 
(|2009|) . who include all of our sample objects except for 
IRAS 13349H-2438, and found similar results. Objects with 
the smallest uncertainties on i?Fe occupy a narrow regio n 
of RpJRg around Ri 20,000km s“^ (cf. iShu et al] 1201 lb . 
There appears to be an intriguing hint that sources with high 
Lboi/Lpm (above Ri0.06) show smaller detected values of 
RpJRg than sources at low Eddington ratios. But this is not 
a signihcant trend in the current sample at least, and for the 
comparison between RpJRg and LBoi/TEdd we find p = - 


0.10^042 (90%). The second panel in the figure illustrates 
that dust is present on scales similar to RpJRg in more than 
half the sample. 

4. DISCUSSION 

4.1. Why does i?dust serve as an outer envelope to i?Fe ? 

Our key result is that the dust sublimation radius forms an 
outer envelope to the Ee Ka emission zone in Type 1 AGN, 
with uncertainties allowing i?Fe at most a factor of a few times 
7?dust for our sample. This is not an obvious prediction given 
that the gas distribution is likely to extend continuously on 
scales both larger and smaller than the innermost dust radius, 
and that tori may be clumpy with low optical depth ‘holes’ 
allowing radiation to penetrate well into its body. Our work is 
thus constraining for AGN torus models, as follows. 

Eor axisymmetric gas distributions, the emitted Ka line in¬ 
tensity at any radius r depends upon th e local solid angle 
^Kn( r) of gas illuminated by the AGN (iKrolik & KallmanI 
Il987h . The absence of sources with i?Fe>>7Jdust 
then argues against geometries in which ^Kair >> 
7?dust) > f^Ka(.Rdust), i-C- where the line emitting area in¬ 
creases strongly with radius beyond i?dust- Line emission will 
not cease exactly at i?dust, which would be an unphysical sce¬ 
nario. Eig.[T] shows that the i?Fe values of up to a factor of a 
few times i?dust are allowed within the uncertainties. How 
do typical torus models match up to these constraints? In 
the geometry adopted by the widely-used MYTORUS model 
(iMurphv & YaaoobI l2009h . the torus is assumed to have a 
donut shape with a circular cross-section and a half opening 
angle of 60°. Radiation from the nucleus that directly illumi¬ 
nates the torus would impinge upon the donut surface out to 
distances of the surface tangent line. In the MYTORUS geom¬ 
etry, this distance is 1.7 i?dust where i?dust is the inner edge 
of the torus. This is consistent within the constraints from 
Eig. 1. On the other hand, tori with strong flaring in height as 
a function of radius would present large illuminated surface 
solid angles at large nuclear distances, and would not easily 
satisfy the above constraints on i?Fe relative to i?dust- 

Radiation would also penetrate below the torus surface, 
but if the integrated column density through the toru^ is 
Compton-thick as is generally thought to be the case, then 
absorption and Compton-scattering would quickly deplete 
line photons in the body of the torus, so Hkq refers to the 
solid angle of the reflecting torus surface down to an effec¬ 
tive optical depth tko ^ L A natural prediction of this sce¬ 
nario is that sources in which the ob scuring tori are likely 
to be Comp ton-thin (e.g. NG C 2110, iMarinucci et ani2015l : 
NGC 7213. lUrsini et alJl2()T^ need not follow the trend of 
Fig.IU i.e. their EeKa emission need not be restricted to lie 
around i?dust- 

We also emphasize that our results do not exclude the pres¬ 
ence of gas on extended scales altogether. If more distant gas 
is not excited by AGN radiation, it will not fluoresce and will 
remain invisible. Eurthermore, studies based upon line width 
measurements will selectively target the strongest emission 
regions present on the innermost (fastest) scales. More dis¬ 
tant (and fainter) emitting components could become appar¬ 
ent if our direct line-of-sight to the inner torus were obscured 
by material optically-thick to the line photons, as may be ex¬ 
pected in Type 2 and Compton-thick AGN (cf. the detection 
of extended emission in NGC 4945 by IMarinucci et alJl2012l 

* It should be noted that this is different from the line-of-sight column 
density, which would vary with inclination angle. 
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^dust (pc) 


Fig. 1.— NIR hot dust sublimation radii (i?dust) vs. Fe Ka line region radii (Rfs) shown in black. Circles represent .Rdust from reverberation lags and squares 
are interferometric measurements. The unfilled violet symbols are the BLR radii (.Rblr) based upon optical reverberation mapping. Dotted lines connect these 
various measurements for any one source. The gray shaded zone represents > i?dust’ scales commensurate with the body of the torus and beyond. 


and in NGC 1068 bv iBauer et Monte Carlo radia¬ 

tive transfer simulations can place detailed constraints on the 
radial gas geometry in the torus. 

There are several potential complexities related to modeling 
of the Fe Ka line and sample selection which must be inves¬ 
tigated to test the robustness of the abov e co nstraints. These 
caveats are discussed at length in Section l44l 

4.2. Fe Ka emission radii much smaller than the torus 

Conversely, the scatter in i?Fe to scales much smaller than 
f?dust argues for Fe K a emission from extended regions. This 
result is not new (cf. lYagoob & Padmanabhanll2004t iNandral 
[2005 lShuetal.1 120101) . but is now inferred from the direct 
size comparison of Fig.[T] The absence of any obvious lu¬ 
minosity scaling in Fig.|2] supports this scenario, and Fig.|3] 
shows that any Eddington rate-driven physical components 
(e.g. accretion disk outflows) do not control the origin of 
FeKa, at least over the range of TBoi/f^Edd that we sample. 
The outer accretion disk and BLR clouds have been invoked 
as Fe Ka emitters in several objects. For example, rapidly 
variable narrow FeKa line from small scal es have been 
found in Mrk509 (iPonti et al.ll2013h Mrk8 41 (iPetrucci et al.l 
200l and NGC 731 4 dVagoob et akllMll . For NGC 7213, 
Bianchi et al.l (l2008ll found that resolved FeKa line width 
matches the opt ical Ha line wid th and argued for a BLR ori¬ 
gin in this case. IShu et all (1201 Ol) found that a signihcant frac¬ 
tion of their sample showed Fe Ka line widths consistent with, 
or broader than, the widths of typical optical and infrared 


photoionized lines, concluding that the Fe Ka emission zone 
location (relative to that of the BLR) genuinely varies from 
source to source. A BLR origin is also consistent for several 
sources in our sample. 

Fig. m shows that in many of the sources with the best con- 
sti-aints on i?Fe (e.g. NGC 4151, NGC 5548, NGC 3783), i?Fe 
is well matched to i?dust- However, this does not imply that 
low grating data S/N is the cause of mismatching i?Fe and 
7?dust values in other objects. It is certainly true that the 
largest FeKa error bars in Figs. [T]and|2] must be unphysical 
and are likely to be a result of the line fits modeling the under- 
lyi ng continuum (se e, for example, the discussion on Fairall 9 
by IShu et al.l 120101) . In addition, some of the sources with 
ill-constrained iipe also have relatively short ^exposure times 
(fexp) - e g- both Fairall 9 and NGC 4593 were observed for 
80 ks with the HEG as compared to a median fexp = 150ks 
for the sample and a maximum of 890 ks in the case of 
NGC 3 783. But as emphasized by I Yaqoob & Padmanabhanl 
(l2004ft . the uncertainties are not simply a function of S/N, 
with the grating data showing that the peak Ee Ka energy of 
6.4keV is not dominated by a narrow core in these sources^ 
Moreover, if we use strict er conhdence in tervals on ufwhm 
( e.g. the 90% range from IShu et alJl20Tol) . the sources with 
i?Fe significantly less than i?dust still remain discrepant in all 

^ Several of these including Fai rall 9 and NGC 4593 be l ong t o the 
‘Group 2’ classification of objects in lYauoob & Padmanabhanl 4200411 with 
weaker or absent naiTow cores. 
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Fig. 2.— Optical AG N luminosity (J/ 5500 ) vs. flpe (fil led) and /Jblr (unfilled). The dotted line is the torus size-lu minosity relation Eg. 1 from 
IKishimoto et al.l <201111 originally from ISuganuma et alJ2006l) and the dashed line is the comesponding relation for the BLR from IBentz et alj <20091) . 


cases. The large uncertainties then appear to reflect an ori¬ 
gin from multiple regions in the AGN environment including 
(fast) clouds present on many scales. Longer grating observa¬ 
tions of the sources with large uncertainties on zifwhm will 
help to confirm the absence of dominant narrow cores, or to 
disentangle fainter core components. 

4.3. Comparison to other recent works 

In two other recent wor ks, iJiang et aH (1201 ll) and 
iMinezaki & Matsushita! (120151) have investigated the use of 
the F eKg line for me asurement of black hole masses in 
AGN. Ijiang et al.l(l201 Ih assumed a torus origin for the Fe Ka 
line with radii derived from infrared reverberation and found 
a consistent trend between the masses predicted assuming 
isotropic motion of the FeKa clouds with the masses esti¬ 
mated from optical reverberation mapping. However, a scal¬ 
ing factor was required for the two sets of masses to match 
each other - the scaling being equivalent to an average line 
core ufwhm that is 2 . 61 q 4 times broader than predicted from 
in itial assumption of a torus origi n. 

IMinezaki & Matsushital (1201 5l) studied a restricted AGN 
sample with ‘best constraints’ on the FWHM of the FeKa 
line, and found that generally the Fe Ka emission zone is lo¬ 
cated between the reverberation radii of the broad H/3 emis¬ 
sion line and the dust torus emission, the latter being esti¬ 
mated statistically (i.e. using a fixed ratio of i?dust to i?BLR 
based on AT-band reverberation measurements). The cri- 
teria that constitute ‘be s t cons traints’ are not quantified by 
IMinezaki & Matsushital (1201 5l) but their selection excludes 
most objects with the broadest Fe Ka lines, as well as sources 
with unresolved FeKa lines. The former category of sources 
includes objects in which the Fe Ka emission zone is consis¬ 
tent with, or sign ifican tly smaller than the BLR (as we dis¬ 
cussed in Section lA^ . while the unresolved sources could 


potentially show F eKa emission fro m outside the torus re¬ 
verberation radius. Lliang et ^ (1201 ll) . on the other hand, do 
include two sources with upper limits on vfwhmO and note 
the possibility that the narrow line core originates from radii 
smaller than the infrared radiation. 

Our main result was derived entirely independent of the 
above works. While being consistent, a key distinction is that 
we go beyond average comparisons on a statistical basis. Im¬ 
portantly, from our object-by-object comparison, we make a 
prediction that i?dust forms an approximate outer envelope 
for normal Type 1 AGN. This prediction can be easily falsi¬ 
fied if i?Fe were observed to be much larger than i?dust in fu¬ 
ture observations (see also discussion in next section). In this 
sense, it is much more powerful than inferences based upon 
statistical trends. This holds d espite the f act tha t we have not 
restricted the base sample of IShu et al.l (120101) according to 
FWHM constraints, and have included upper limits as well. 
We have then explored the consequences of such a postula¬ 
tion for the structure of AGN tori (Section 4.1) - aspects not 
mentioned in these other wo rks whose focus is Mrh mea- 
surement. We note that since IMinezaki & Matsushital (120151) 
do not include sources in which the peak Fe Ka energy is not 
dominated by a narrow coreQ they do not probe the smallest 
scales that we discuss in Section 4.2 above. 

Next, we will address various potential caveats which could 
affect our main result. 

4.4. Potential caveats 

Several potential caveats need to be kept in mind regard¬ 
ing current analyses of ^Fe Ka line profiles, and these are dis¬ 
cussed in detail below. 

^ They treat Mrlc 509 as an upper limit on ttpwHM. The resultant lower 
limit on ijpg is entirely consistent with our trend between Rdust and Rpe ■ 

^ Except for their one “outlier” NGC7469. 
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Fig. 3.— (Top) Eddington ratios (TBol/TEdd) vs. ijpe in units of Rg. (Bottom) Tsol/TEdd vs. iJdust in units of Rg. Symbols are as in Fig.[T] 


4.4.1. Combining multiple HEG datasets 

For several sources (e.g. NGC3783, NGC4051) multiple 
HEG datasets have been combined in order to maximize the 
signabnoise and place the best constraints on the line width. 
This methodology ignores any potential change in the line 
emission zone between obs ervations. From the individual 
measurements presented by IShu et al] (1201 Ol) . there appears 
to be some variation in FWHM(Fe Ka) for 2 of the 5 sources 
with multiple observations spaced by a few months to years 


(these 2 sources are NGC3516 and NGC5548). But this is 
true only at the 68 % confidence level, and the 90 % con¬ 
fidence intervals show no such variation. In addition, the 
mean FWHM(Fe Ka) measured per sour ce (when combin ing 
multiple observations) for the sample of [ShueFal] <1201 Oh is 
entirely consistent with the mean FWHM(Fe Ka) measured 
per observation. These facts suggest that combining multi¬ 
ple observations does not strongly affect our results. How¬ 
ever, note that for two sources from our sample (3C273 and 
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IRAS 13349+2438), the Fe Ka line is detected in only one ob¬ 
servation and it is unknown whether this is due to intrinsic and 
strong variability of the line. 

The effective FWHM in a combined dataset can also be ar¬ 
tificially broa dened if the cen troid energy changes between 
observations. IShu et all (1201 Oh find a sharp peak in the rest- 
frame centroid energies for their sample at ^ 6.400 keV and 
quote mean sample line energies centered on the neutral line 
to within -12 eV and +3 eV (in comparison to the HEG spec¬ 
tral resolution of ^40 eV), with the per source analysis fully 
consistent with the per observation one. Only in 2 cases out 
of 32 (3C 273 and 4C 74.26) did the centroid deviate strongly 
from the expected neutral line energy, although they note that 
these cases are marginal. 

We also note that the above centroid energy range is equiv¬ 
alent to systemic velocity shifts of r; 140 kms“^ (blueshift) 
up to « 560 km s“^ (redshifts) for the fluorescing material. In 
contrast, the median escape velocity at i?dust for our sources 
is ^2400kms“^ and the median line width for sources with 
an FWHM detection is 3400km s“^, arguing against any bias 
related to systemic outflows or inflows. 

4 . 4 . 2 . The limited HEG spectral resolution 

What is the impact of the limiting spectral resolu¬ 
tion of HEG for our results? The velocity resolution 
of ufwhm = 1860km s“^ corresponds to i?Fe = 0.08pc for 
ATbh = 5x 1O^M0 (the median Mbh for our sample). 
Whereas much narrower lines can be deconvolved given 
enough S/N in the data (e.g. see the case of NGC 3783 which 
has the best S/N and a resolved line even at 99 % confidence), 
this may not be possible if the S/N is weaker. So are we bi¬ 
ased to objects with broader lines, hence creating an artificial 
limit to i?Fe? We argue that this is not the case, because for 
any given line flux, it is always easier to detect the line if it is 
unresolved. If the line is significantly broader than the spec¬ 
tral resolution, its peak decreases and the flux is spread over 
a range of energy bins, in which case detection naturally be¬ 
comes more difficult. Despite this fact, most sources in our 
sample already show resolved lines with the HEG. And we 
emphasize that we are not excluding objects with unresolved 
lines which would show only upper limits to ufwhm- To¬ 
gether, these argue against an artificial limit being created by 
the HEG spectral resolution in our analysis. 

4 . 4 . 3 . Biases due to multiple unmodeled line components 

Another caveat is related to modeling the 6.4 keV feature 
as a single Gaussian and ignoring various intrinsic contri¬ 
butions to the feature such as the Compton shoulder which 
results from down-scattering of line photons off circumnu- 
clear gas. However, the flux in this component is expected 
to be only ^10-20% of that in the narrow line core for re¬ 
flection off Compton-thick material (lMatdl20d^ . so a single 
Gaussian fit would only be slightly broader than in the ab¬ 
sence of a Compton shoulder. We quantified this effect by 
simulating HEG Ee complex spectra for a range of column 
densities upto Nu = 10^^ c m~^ using the MYTORUSL model 
(iMurphv & YaqoobI l2009tl . which self-consistently models 
fluorescence and Compton scattering off a toriodal reflector. 
Another set of simulations using the GAUSS model alone for 
a single line were carried out for comparison. The intrinsic 
Ee Ka line width was set to cr = 1 eV, well below the HEG 
resolution. We used the standard first order HEG canned re¬ 


sponse matrices and effective area flle^ for the simulations 
and covered a broad range in S/N similar to those in the real 
data. All spectra were then fit with single Gaussian lines. We 
found that the effective Gaussian line widths for the simulated 
MYTORUSL spectra were at most 18% broader than for the 
single GAUSS models. Correcting for this would increase 
by at most 39% or 0.14dex - at most a minor effect as can 
be seen in Eig.[T| An ind ependent constraint comes from the 
fact that lShu et al.l (1201 Oh found a sharp centroid peak for their 
sample of close to the expect mean rest-frame line energy at 
6.400keV. If the Compton shoulder were strongly skewing 
the line fits, the centroid would be expected to lie at lower 
energies, but this not the case. 

Similarly, blended lines could produce artificial broadening 
in single Gaussian fits. Such lines may arise from multiple 
ionization stages of Ee, and we here quantify their influence 
on the measurement of the effective EWHM, and hence on the 
derivation of i?Fe- Eirstly, it is worth emphasizing that the ob¬ 
served line centroid energy distribution peaking at r: 6.40 keV 
for the full HEG sample argues strongly for low ly ionized 
states (typically below Ee XVII), as discussed by IShu et al] 
(1201 Oh. Upon exam ining the centroid energies (E) reported 
bv IShu et al.l (1201 Oh for our sample of 13 AGN, we And a me¬ 
dian i? = 6.403keV (iO.OlOkeV mean standard deviation). 
The upper end of this range (ii/ = 6.413ke V) corresponds 
closely to an ioniza t ion stage of Ee XIV (c f. iLiu et al.ll2010t 
iPalmeri et al.l 120031 : iMendoza et alJ 1200^ . A further con¬ 
straint on the strongest ionization stage comes from the distri¬ 
bution of E values for the resolved lines in the longest HEG 
exposures. There are 5 objects with exposure time of more 
thanR! 150ks: NGC3516,NGC3783,NGC4151,NGC5548 
and IRAS 13349+2438. All except IRAS 13349+243^ show 
£^<6.415keV at the conservative 90% confidence level, 
again consistent with ionization stages lower than Eexiv- 
XV. Based upon this upper energy range, we assess a worst 
case scenario where all ionization stages from Ee I to Ee XVIV 
contribute maximally to artificially broadening the line. All 
stages are taken to possess the same component line strength 
and the same intrinsic EWHM. The centroid energies for 
these stages ran ge over i? r; 6.394-6.415keV (cf. Eig.2 of 
ILiu et al.ll20l0h . or a full range of A£' = 21eV. This corre¬ 
sponds to a velocity broadening of 980km s“^. This velocity 
should be subtracted in quadrature from the measured EWHM 
values when modeling the 6.4keV feature as a single Gaus¬ 
sian. The largest effect will be on the source with the smallest 
EWHM for which the line is resolved, which is NGC 3783 
(EWHM= 1750kms“^; see Table 1), and in this case, the in¬ 
trinsic EWHM of the individual ionization stages is only 20 % 
smaller, at 1450kms“^. The corresponding i?Fe increases by 
45 % and is still fully consistent with i?dust- The effect of line 
blending for the other sources in our sample is much less. 

4 . 4 . 4 . Sample representativeness 

The parent sample of IShu et al.l (120101) is certainly not a 
complete sample in a physical sense that is ideally suited for 
statistical studies. Despite being the largest sample with high 
quality Ee Ka measurements to date, there is some unavoid¬ 
able bias of HEG observations towards sources with a previ¬ 
ously known presence of Ee Ka (either broad or narrow) and 

^ http : //cxc . harvard. edu/caldb/prop^jlan/grating 

® This source is already noted bv IShu et alj 420101) as one of the few anoma¬ 
lous cases with potential line variability between observations. The data for 
this source also yield only very weak consteaints on the overall line properties 
including flux and FWHM (see Fig. 1). 
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to objects that were observ ed with other aims in mind (see also 
discussion on th is point in lYaaoob & Padmanabhanl2004l and 
iShu et al.ll2()I^ . 

On the other hand, we note that there is no bias in terms 
of physical source properties or classification. The sample 
includes optical Type 1, 1.5 and 2 AGN and also some Nar¬ 
row Line Seyfert Is, all below z = 0.3. However, the sam¬ 
ple excludes obscured sources with line-of-sight gas column 
densities A^h > 10^^ cm“^ in order to avoid spectral modeling 
complexities. While a signal-to-noise requirement is imposed 
so that the t otal count s in the full HEG bandpass must be more 
than 1500, iShu et al] (120101) state that relaxing this criterion 
would have admitted only 2 more objects, so this does not ap¬ 
pear to be a major restriction unless Fe Ka properties evolve 
very rapidly with redshift (i.e. to fainter source fluxes) - this 
is an unlikely prospect in our op inion. _ 

Of the 36 sources studied by IShu et al] (|20T3), the statis¬ 
tical quality of the data was good enough for FWHM mea¬ 
surements in 27 objects. In the remaining 9 objects, the sig- 
naknoise of the data around the Fe line region is so poor 
that the line itself is not detected at more than 95 % confi¬ 
dence and the Gaussian model fit can become unphysically 
broad as it begins to mo del the continuum itself (see Sec¬ 
tion 3 of IShu et al.ll201(i) . This is the case for NGC526a, 
Mrk705, NGC3227, Mrk766, PDS456, IRAS 18325-5926, 
NGC7314, Ark564 and MR2251-178. Examining the rea¬ 
son for the poor data quality of these 9 sources, we note that 
their mean HEG exposure time is only 85 ks, as opposed to 
192 ks for the 27 objects with FWHM measurements. In 
addition, the two sources with the longest exposure times 
(PDS456 and MR 2251-178 with exposures of 145.2ks and 
148.7 ks, respectively) are both known to be relatively dis¬ 
tant, obscured quasars showing highly complex continua and 
evidence of sign ificant absorption lines resulting from st rong 
outflows (e.g. iGibson et alJl200-'^ iHagino et al.ll20L5l) . In 
other words, the absence of a significant Fe Ka emission line 
in these 9 sources is most likely a combination of relatively 
short exposure times resulting in low signaknoise, and com¬ 
plex spectra in which a uniform analysis does not immediately 
resolve a narrow line core. 

The main restriction imposed by the infrared selection is 
one of target visibility - most (though not all) of the dust re¬ 
verberation lags have been carried out from the MAGNUM 
telescope in Hawaii. There may be an implicit preference for 
observing bright sources and those sources with known (sig¬ 
nificant) levels of optical variability, but it is not clear how this 
would bias our sample. Interferometric size measurements 
have been carried out mostly from Keck also in Hawaii, but 
the sample of AGN where i?dust,intf has been measured is 
currently limited. So there is no obvious physical bias in the 
source selection, albeit our sample of only 13 AGN from the 
parent sample remains small. 

4.4.5. Alternate black hole mass estimates 

Finally, one may question the robustness and suitabilit y of 
our adopted Mbh measurements. As discussed in Section im 
certain assumptions about the geometric correction (/-factor) 
are implicit in virial Mbh estimates based upon reverberation 
mapping, and we have assumed that similar correction factors 
apply for the Fe Ka emission zone. 

The influence of making such an assumption can be judged 
by using alternate black hole mass measurements, which are 
entirely independent of reverberation mapping. In the Ap¬ 
pendix, we have compiled such measurements in Table|2] We 


present a figure showing the resultant comparison of i?Fe with 
^dust in this case, and demonstrate that our main inference of 
^dust forming an envelope to i?Fe remains unchanged. We 
refer the reader to the Appendix for this test. 

To summarize, we have discussed and quantified some 
caveats to our work, including those from combining multiple 
datasets, limited instrumental resolution and biases related to 
simplistic (single Gaussian) modeling of a likely complex line 
feature. None of these issues appears to cause a bias strong 
enough to overwhelm our main result of i?dust acting as an 
outer envelope to i?Fe- We point out current shortfalls in terms 
of sample completeness, although the apparent heterogeneous 
nature of sample selection argues against any strong selection 
bias. We also explore the use of alternate black hole masses 
and their influence on our main result, which do not change 
our main result. 

Although we cannot rule out the presence of a combination 
of these various biases, it would be surprising if they were to 
conspire to produce a limit of i?Fe around i?dust - which is a 
completely independent quantity. However, we acknowledge 
the current limited sample size and the paucity of line variabil¬ 
ity measurements. More complete sampling of AGN in both 
X-rays and in the infrared are clearly important for drawing 
robust conclusions, and we discuss future possibilities in the 
next section. 

4.5. Future perspectives 

This held is expected to lea p forward with the im minent 
launch of the Astro-H mission (iTakahashi et akll2014ft . With 
a spectral resolution of ~ 4-7 eV, the Soft X-ray Spectrome¬ 
ter (SXS) employing calorimetric photon detection will pro¬ 
vide an improvement in spectral resolution of Ri 6-10 as com¬ 
pared to present grating instruments. Since Rpe has an inverse 
quadratic scaling with line velocity, SXS will extend the ra¬ 
dius out to which the Ka emission can be localized by a factor 
of ^ 40 in any individual object. This assumes a single Gaus¬ 
sian feature. Individual ionization stages will not all be neces¬ 
sarily separable by the SXS, though some deconvolution will 
be possible. In terms of the overall sensitivity to line parame¬ 
ter measurement, this can be quantified in terms of a ‘figure of 
merit’, which combines the effective area and the spectral res¬ 
olution. Relative to HEG, the overall SXS sensitivity to line 
detection is expected to be improved by a factor of r; 7, ex¬ 
tending our reach to correspondingly faint (and more distant) 
systems. However, the improvement in the corresponding fig¬ 
ure of merit for detection of line broadening (which is key 
to our work) is about two orders of magnitude relative to the 

hegQ _ 

In the next decade, the Athena mission (iNandra et alJ201^ 
will have a much larger collecting area than Astro-H, pushing 
such studies out to high redshift. Reverberation of the nar¬ 
row FeKa core for more objects would also be an indepen¬ 
dent constraint on Rpe. Variable narrow lin es in som e cases 
have already been mentioned. In addition, iLiu et alJ (120101) 
found Fe line reverberation on a timescale of ^ 20^0 day in 
NGC5548. This would push Rpe a factor of ^ 1.5-3 lower 
than our estimate in Table 1, but still consistent within the un¬ 
certainty on Rpe. 

And in a similar vein, enlarged NIR coverage will greatly 
help to better All the Rpe vs. i?dust parameter space. There 
are currently a total of 31 AGN with measurements of i?dust- 

^ http://astro-h.isas.jaxa.jp/ahqr.pdf 
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This includes sources from the references cited in Table 1 as 
well as a few other publish ed and unpublished sources (e.g. 
GQCom; ISitko et al]ll993h . We have included 13 of these 
(42%) in our X-ray cross-matched sample. Although Keck 
Interferometer is no longer available, more size measurements 
will be available soon (Kishimoto et al. 2015 submitted) mak¬ 
ing use of the three-telescope beam combiner AMBER at the 
VLTI. This will boost the present sample and address pro¬ 
jection effects by invoking closure phase data. And although 
the MAGNUM telescope is now decommissioned, there are 
many other worl d-wide efforts ongoing t o obtain AGN NIR 
time lags (e.g. iPozo Nunez et alJl2015h . LSST could also 
serve as a dust reverbera tion machine providing thousands of 
time lag measurements (lHonidl2014t) . Finally, in the distant 
future, we can expect to directly resolve the torus in many 
more sources with larger interferomete rs such as the Planet 
Formation Imager (iMonnier et al.ll2M4h . 
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P.G., a Marie Curie International Incoming Fellowship within 
the EC Framework (PIIF-GA-2013-623804) for S.F.H. 
and JSPS (26887044) for M.K. We thank C. Done, Y. Ueda, 
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tonucci for comments on an initial draft which greatly im¬ 
proved the discussion of the robustness of our results. We 
also acknowledge useful comments from the anonymous ref¬ 
eree on the first submission and the revised draft. 
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^dust (pc) 

Fig. 4.— The comparison of Rpe with i?dust,rev in Fig. 1, but this time with the computation of R^e based upon black hole masses obtained from 
reverberation-independent methods. 


APPENDIX 


BLACK HOLE MASSES INDEPENDENT OF REVERBERATION MAPPING 

Our adopted values of Mbh in Table 1 are mostly virial estimates based upon reverberation mapping of the BLR, which 
implicitly assume a certain correction for geometric projection effects. The size comparisons in Fig. 1 effectively presuppose that 
the same projection effects apply for the Fe Ka emissio n reg ions and the BLR, which may, or may not, be the case. Although the 
correction factors are expected to be small (see Section lOT . it is worth investigating how strongly the above assumption affects 
our results. In order to assess this, we compile black hole masses which are independent of reverberation mapping. Any distinct 
projection effect for the Fe Ka and H/3 emission regions will then be retained in the emitter size comparisons, and the resultant 
change in the i?Fe distribution (relative to Fig. 1) will give an estimate of the maximal impact that the unknown projection effects 
can have. 

Table|2]presents our compilation of alternate Mbh values0 The se are mostly b a sed up on host galaxy stellar velocity dispersion 
(cr,) measurements and the relation between Mbh and fromlGiiltekin et al.l (l2009h. W here this was not possible, we used 
scaling relations between X-ray variability amplitude (IZhou et 120101: iMcHardvl 12013h and Mbh, or the empirical relation 
between optical continuum luminosity and the BLR radius dVestergaard |200^"combined with the FWHM of the H/3 emission 
line. There is broad agreement between the alternate masses tabulated here and those listed in Table 1, with the alternate masses 
(mostly based upon Mbh-ct*) being larger than the virial estimates by a median factor of 1.8 but with a significant spread of 0.7 
dex (standard deviation) between the two sets of masses. 

The resultant values of i?Fe are also listed in the table, and are plotted in Fig.|4]which shows that although the location of the 
FeKa emission radii relative to i?dust (and thus also relative to TJblr) do change in some individual objects as compared to 
Fig. 1 (e.g. i?Fe(NGC 3783) is now smaller than i?dust, and i?Fe(Mrk590) is now consistent with i?dust), our main inference of 
7?dust serving as an approximate envelope to i?Fe still holds. IC4329A is t he o nly source in which the limit to i?Fe apparently 
sits very close to (but still consistent with) i?dust- As we discussed in Section l4n Fe Ka emission will not cease exactly at i?dust, 
and i?Fe values of up to a few times i?dust are allowed by Fig. 1 and also by typical compact torus models. Astro-H observations 
of this source will help to pinpoint its location on the i?Fe-7?dust plane. 
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TABLE 2 

Alternate black hole mass estimates 


Source 

Method(MBH) 

Measurement 

Ref. 

log Mbh 

AlogMBH 

Rpe 





Mq 

Mq 

X 0.1 pc 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Fairall 9 

Mbh-ct* 

cr+=228± 18 km s”'*' 

Oliva et al. ^99) 

8.36±0.46 

+0.05±0.47 

U.UH-W_Q Q4Q 

Mrk590 


fT*=189±6kms“^ 

^jelsor^^^2004) 

8.02±0.44 

-0.82±0.45 

0 '214+0-890 
u.ji't — Q 298 

NGC3516 

Mbh-<7* 

fT*=181±5kms“^ 

Nelson et al. 12004) 

7.94±0.44 

-0.31±0.47 


NGC3783 

Mbh-<7* 

o-,=95±10kms“i 

^nkei^^^20M) 

6.75±0.48 

+0.73±0.49 


NGC4051 

Mbh-<7* 

fT*=86±3kms“^ 

Nelson et al. ^2004) 

6.57±0.44 

-0.28±0.47 

>0.061 

NGC4151 

Mbh-<7* 

fT*=116±3kms“^ 

^nkei^^^20M) 

7.12±0.44 

+0.57±0.48 

0 148+0-288 
u.ih-o_q 098 

3C273 

X-ray variance 

log(NVA) = 0.20±0.08 

^cHa^^20^) 

8.86±0.08 

+0.96±0.12 

i.mtllll 

NGC4593 

Mbh-<7* 

fT*=135±6kms“^ 

Nelson et al. ^2004) 

7.40±0.45 

-0.67±0.88 

0.2871°;™? 

IRAS 13349+2438 

Lbioq-Rbbr 

srg s ^ 

Wans et al. (2004); Grune et al. (2004) 

7.74±0.50 

+ 1.26±0.66 

0 119+0-474 

O.llV_0 119 

IC 4329A 


(T.=225 ±9kms“^ 

^liv^^^^^99) 

8.24±0.44 

-1.25±0.90 

> 2.911 

NGC5548 
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Woo et al. (2010) 
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1 0981;^;"? 
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Mbh-cT:^ 

cr*=131±5kms“^ 

Nelson et al. (2004) 

7.34±0.45 

-0.25±0.45 

0.053i“;J™ 


Colurrm (2) denotes the method used to determine Mnn. Three methods are used: the relation from lGiiltekin et al.ll20Q9l) : the X-ray variance scaling with Mt^vi either 

from lMcHardvl I2Q131 who characterizes this in terms of the Normalized Variability Amplitude or NVA) or from IZhou et alj <20101 who provide measurements of the rms variability 
(7^ms)l empirical relation of L 5100 with -Rblr from IVestergaarJ l2QQ2h . Column (3) lists the relevant observable used in the method. Column (4) gives the reference for the 
measurement of the observable. Column (5) lists the resultant ATbh value, and Column (6) the difference in log Mbh with respect to the value quoted and used in the main paper 
(Table 1). Finally, Column (7) lists the corresponding value of R^e using the alternate Mbh estimate. For IRAS 13349+2438, we assume a large error of 0.5 dex in Mbh as in 
IDe Marco et aljf2Q13I) . Note that for this source Up^j^j^=2800 ib 180kms“^ is reported in IGrupe et al J f2QQ^ . 
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